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Abstract 

We discuss the theory of knots, and describe how knot invariants arise naturally in 
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jy-^ gravitational physics. The focus of this review is to delineate the relationship between 

knot theory and the loop representation of non-perturbative canonical quantum general 
relativity (loop quantum gravity). This leads naturally to a discussion of the Kodama 
wavefunction, a state which is conjectured to be the ground state of the gravitational field 
with positive cosmological constant. This review can serve as a self-contained introduction 
to loop quantum gravity and related areas. Our intent is to make the paper accessible to 
| a wider audience that may include topologists, knot-theorists, and other persons innocent 

of the physical background to this approach to quantum gravity. 
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1 Introduction 

"Physical laws should have mathematical beauty." ~P AM Dirac, Moscow 1955. 

The oldest known invariant of knots is the Gauss linking number, which was discovered 
by Johannes Gauss [1,2] in 1833. Physically this is the work done in transporting a magnetic 
monopole around a closed loop L\ through a magnetic field that is produced by a current 
moving through a closed loop of wire Li- The electric current and magnetic field have fixed 
directions, so this configuration may be thought of as an oriented link with two components. 
Topologically the invariant measures the degree to which the two loops are linked. Furthermore, 
the invariant is well defined so long as the two loops do not touch or intersect. Gauss calculated 
the work done on the monopole to be the double line integral 
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The invariant is then the work divided by the current I and magnetic charge m of the monopolc. 
A generalization of this invariant was given by Edward Witten [3] in 1989 as the expectation 
value of certain gauge-invariant functionals, showing that topological invariants of links are 
fundamental in gauge theories. 

The analytical study of knots began in the 1880's, after James Maxwell [4] had established 
the equations of electrodynamics on a firm mathematical foundation. At that time quantum 
mechanics had not yet been postulated, and the remarkable stability of atoms was a mystery. It 
seemed quite natural therefore to explain these structures purely in terms of Maxwell's theory. 
William Thomson (Lord Kelvin) [5] put forward the hypothesis that atoms were knotted flux 
tubes of the aether. His motivation really came from the strong analogy between the vacuum 
equations for the electric field and the fluid flow equations of hydrodynamics, but this also gave 
a satisfactory explanation for the stability and large variety of atoms. This theory was taken 
seriously for a while and attracted a modest following. In particular, Peter Tait [6] began his 
classification of knots as projections onto a two-dimensional plane. Along the way he noticed 
some properties which led him to what we now call the Tait conjectures. Most of these were 
not proved until the 1980's after Vaughan Jones [7] discovered a polynomial invariant which 
now bares his name. Indeed this discovery marked a revival of knot theory long after the 
inevitable failure of the vortex model when the principle of relativity by Albert Einstein [8] did 
away with the aether, and with it the possibility of vortices in the aether. In the meantime, 
(up until 1984) knot theory had developed into a vigorous branch of low dimensional topology, 
using all the techniques of modern geometric and algebraic topology. The advent of the Jones 
polynomial and its generalizations brought new techniques and startling new relationships with 
physics into this already active field. After Jones' discovery, many papers began to appear 
in which new link invariants were defined. A definition of the Jones polynomial was given 
by Louis Kauffman [9] as a state summation in analogy with partition functions in certain 
low dimensional systems in statistical mechanics (such as the Potts model). This model is 
called the bracket polynomial or the Kauffman bracket. The Jones polynomial is the product 
of the bracket polynomial and the self-linking number (the writhe) of a given link. Seminal 
work of Witten [3] then showed that link invariants could be constructed via gauge theory. 
Witten accomplished this relationship by calculating the expectation values of Wilson loops 
using facts from two-dimensional conformal field theory. Simultaneously and independent of 
this result, a representation for four-dimensional canonical quantum gravity was proposed by 
Carlo Rovelli and Lee Smolin [10, 11], which exploits the topological invariance of physical states 
in the so-called loop space, and are otherwise invariant under spacetime diffeomorphisms. This 
generalized earlier work of Ted Jacobson and Smolin [12], who showed that loop solutions solved 
all the constraints of the Ashtekar phase space. Rovelli and Smolin were motivated by an idea 
due to Chris Isham [13], that a quantization of general relativity may involve a non-canonical 
Poisson algebra. Remarkably, however, this loop representation is closely related to an old idea 
of Roger Penrose [14], whereby the fundamental structure of spacetime is combinatorial, given 
by the spin networks. Soon after, many exact states satisfying the quantum constraints were 
found by Hideo Kodama [15,16], Viqar Husain [17] and Bcrnd Brugmann, Rodolfo Gambini 
and Jorge Pullin [18]; these were all shown to be related to the Jones polynomial and other 
similar knot invariants. As we will see, none of the ideas presented here are independent of each 
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other. 



The outline of the paper is as follows. We begin §11 by reviewing the foundations of non- 
Abclian gauge theories. We define the holonomy of the connection and the gauge invariant 
Wilson loop. We give a detailed study of the Chern-Simons gauge theory on a three-dimensional 
manifold. In particular, we show how the WZW action is induced as a boundary term under 
a gauge transformation of the Chern-Simons theory, and how the wavefunctions in the Hilbcrt 
space of the Chern-Simons theory are recognized as the conformal blocks that determine the 
iV-point correlation functions of the WZW model. In §111 we move on to knots and knot invari- 
ants. We state the definitions that are needed for the classification of topologically incquivalent 
closed curves embedded in three-dimensional Euclidean space R 3 , and present the topological 
manipulations that are permitted on the projections of knots onto a two-dimensional plane (the 
Reidemeister moves). We then define some of the relevant invariants - the linking number, 
self-linking number (writhe), the Jones polynomial and the bracket polynomial. By defining 
the partition function for an arbitrary manifold as the (unnormalized) expectation value of a 
Wilson loop with the exponential of the Chern-Simons action as the measure, we are able to 
derive the skein relation for the Jones polynomial, showing that the invariant has an intrinsic 
three-dimensional definition. We also derive the Gauss linking number which arises when the 
gauge group is U(l). This completes the background that is neccessary to understand the phys- 
ical states of the quantum gravitational field. In §IV we describe the second-order (metric), 
first-order (tetrad) and self-dual (connection) formulations of general relativity with a non-zero 
cosmological constant. We use the latter with a densitized triad as conjugate momentum to 
derive the phase space which contains seven primary constraints that vanish weakly on the 
constraint surface. We define the fundamental loop observables of the quantum theory and 
describe the spin networks that arise in loop space. We end the formal discussion with the Ko- 
dama state, in particular how it is simultaneously a semiclassical state as well as an exact state, 
how it is equivalent to the bracket polynomial of framed links in loop space and how the SU(2) 
spin networks are q-deformed due to the requirement of a framing of the loop observables. 

2 Chern-Simons Theory and Conformal Field Theory 
2.1 Loop States in Gauge Theory 

Let us begin by reminding the reader of the basic definitions of non-Abelian gauge theories. We 
have a vector potential (connection) one-form A = A^T a dx^ that is valued in a Lie algebra with 
gauge group G and generators T a . This connection determines the gauge covariant derivative 
and curvature by 



The curvature two-form F = {l/2)F^ v T a dx^ A dx v is defined on an arbitrary four-dimensional 
manifold M. . A gauge transformation g e G is a zero- form whose action on the connection and 
curvature is given by 



Fpv = [£>„,£>„] =0„4 



d v A^ + \Af l , A v \. 



(2) 
(3) 



A : A =g~ 1 dg + g~ 1 Ag and F : F = g~ 1 Fg. 



(4) 
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The space of all connections is denoted A, and the set of all gauge transformations is said to 
form the group Q. We write A/Q for the space of gauge equivalent classes of connections, i.e. 
the space of all connections modulo gauge transformations. The action for pure Yang-Mills 
theory [19, 20] on M is 

S YM [9^,A] = l [ Tr(FAF*) = -][ d 4 xJ\det( 9tlv )\Tr(F^F^)- (5) 

Here Tr is a non-degenerate invariant bilinear form on the Lie algebra called the "metric" 
j ab = pTr(T a T b ) for some normalization p and F* is the Hodge dual of F. Variation of the 
action with respect to A yields the equations of motion DF* = 0. 

Now let 7 = -f(s) be a curve that is parametrized by s. The path-ordered product of the 
connection A with components = A^(x(s)) e A is defined such that 

PA^xis!)) ■ ■ ■ A^xisn)) = A AJ (a;(si))---A AJ (a;(s„)), si > ■ ■ ■ > s„. (6) 

This is a permutation of factors ordered so that larger values of Si appear left to right. Now 
define on 7 the path-ordered exponential (called the Wilson line) 



Wfr] = Pexp 



•L 



ig / dx^A^ 
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(7) 



This is the holonomy of the connection A^ on the principal G-bundlc. It is required to transform 
as W[j] — > U(z)W[y]U^ (y) under a general gauge transformation U. If we take 7 to be a closed 
loop then the loop integral transforms as W[j] — > U(y)W[y]U^(y), and taking the trace of this 
gives TrW{-f] -> Tr (U(y)W[i\U*(y)) = TrWfy]. This implies that the relevant quantity that 
is gauge invariant is 



W[j] = TrPexp 



ig (h dx^Af, 



(8) 



called the Wilson loop [19]. If G = U(l) then the path ordering and trace are irrelevant since 
the potentials commute. In this case the Wilson loop gives a measure of the phase that is 
acquired by a charged particle as it moves around the path 7 through an electromagnetic field 
given by = d^A v - d v A ll . 

Let us now introduce some notation for the Wilson loops that we will be needing below. If 
7i represent segments of curves then 7j o 7^ says that the beginning of 7^ is glued onto the end 
of 7j. The segments ji and jj may not be closed, but the composition of the two must form a 
closed loop. The segment A is necessarily an open curve. Also the inverse 7 -1 is the segment 7 
with opposite orientation. With this, the Wilson loops satisfies the following properties [21, 22]: 

W[0] = 1 (9) 
W[ 7 ] = W^oAoA- 1 ] (10) 
^[71072] = W[7 2 o7i]. (11) 

In addition, the Wilson loops also satisfy two properties which hold only for gauge group 
G = SL{2,C). These are 

W[ 7 ] = WIj- 1 } (12) 
W[7i]W[ 72 ] - i(W[ 7 ioAo 72 oA- 1 ] + >V[7ioAo7 2 - 1 oA- 1 ]), (13) 
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The properties (9)-(13) are called the Mandelstam constraints. Diagrammatically for the simple 
case where 71 ~ -Q and 72 ~ -Q (13) says that 

wio ]w[0] = ^ w[ €0] + w [-OQ])- (14) 

2.2 Chern-Simons and Wess-Zumino-Witten Actions 

Let A — AfT a dx l be a connection one-form valued in a Lie algebra with gauge group G and 
generators T a . The Chern-Simons action on an arbitrary three-dimensional manifold A^3 is 
defined as 

S CS [A] = Tr^AAdA+^AAAAA^j 

= A f e vk Tr ( Ai{djAk _ dkA ) + Afe] 



(15) 



where fc is a coupling constant called the level. The equations of motion obtained from variation 
of the action with respect to A are F = dA + A A A = 0, which means that the field strength F 
vanishes. If g e G is a gauge transformation defined by (4), the Chern-Simons action transforms 
as [23,24] 



S CS [A] = S CS [A] - A / Tr Udgg' 1 ) A A 



dM 3 

1 TT[(g- 1 dg)A(g- 1 dg)A(g- 1 dg)]. (16) 
M 3 
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If AI3 is compact and without boundary then the boundary term above vanishes, and Scs is 
gauge-invariant (modulo 2ir) if k is an integer. The last term is an integral multiple of 2ir called 
the winding number of g. 

If M3 has a boundary then a surface term needs to be added to the action to ensure that 
Scs bas proper extrema. This is accomplished by imposing boundary conditions on a field, i.e. 
a particular field is held fixed at dM 3 . For two-dimensional manifolds there is a remarkable 
property that conformal structures are the same as complex structures. See, for example, 
Chapter 21 of Hatfield's book [19]. It is then natural to fix a complex structure on dM^ by 
introducing a complex coordinate z = x 1 +ix 2 with conjugate z = x 1 — ix 2 if the real coordinate 
system is given by (x 1 ,^ 2 ). Choosing the axial gauge A = results in a connection that is 
two-dimensional. The components of this connection transform as 

A 1 = A,gL + A,£L = A, + A,, A a = A,£ + A,£L=i{A,-A,) (17) 

so that in the holomorphic representation (in terms of z and z) the connection is written as 
A = A z dz + Agdz. Fixing the component A z on dMs gives the action 

S$ S [A} = — [ Tr \A AdA+ -A A A A A) + — [ dz A dzTi(A z A z ). (18) 
47r J M 3 V 3 / 4?r J aMs 

Alternatively, fixing the component A z on dMs gives the action 



S® [A] = A f Tr \A AdA+ -A A A A A) — / dz A dzTr(A z A E ). 

47r J Ms V 3 / 4?r J dM3 



(19) 
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When the action (18) or the action (19) are varied, then the boundary terms cancel as can 
be verified. Now, under the gauge transformation g G G (defined by (4)), the action (18) 
transforms as 

SBsW = S§ s [A} + kS® zw lg,A z ] (20) 
S% zw [g,A z ] = -L / T I (g-^d z gg- 1 d z g-2g- 1 d z9 A z ) + -^ f Tr^dg) 3 . (21) 

47r JdM 3 iZ7r J M 3 

where S^zw is the "chiral" Wess-Zumino-Witten action [25,26] on the boundary (see 
below for definition of chirality). This establishes, classically, the connection between Chcrn- 
Simons gauge theory in the bulk manifold and conformal field theory on the boundary. 



2.3 Conformal Field Theory 

We will now discuss some general features of two-dimensional conformal field theory Our 
treatment here is necessarily brief. The main purpose is to introduce the conformal blocks 
which are essential for the quantization of the Chern-Simons theory discussed above. We follow 
the review article by Gaberdiel [27], to which the reader is referred for more details. 

In rf-dimensional Euclidean space R d , with d > 3, the conformal group consists of rotations 
and translations that preserve angles and lengths, the scale transformations given by 

x» -> ilx 11 , G R d , Qe M, (22) 

and the so-called special conformal transformations given by 

G K d , x 2 = x^, a 2 = a^. (23) 



x 2 a^ 



l + 2{x-a) + x 2 a 2 ' 

For d = 2, however, we can introduce the complex coordinates z,z for (a; 1 , a; 2 ) G M 2 such that 
z = x 1 + ix 2 and z — x 1 — ix 2 as we did for above. Here a conformal transformation is 

defined simply by the analytic map z — > f{z) for some (locally) analytic function f(z). The 
conformal group in two dimensions is therefore determined by the set of all analytic maps from 
the plane to itself, which is infinite-dimensional. A subgroup of such mappings form the set of 
Mobius transformations given by 

z — > f(z) = aZ ^j , a,b,c,deC, ad-bc=l. (24) 
To each such transformation can be associated a matrix 

A = ^ ^ , det(A) = 1. (25) 

This group is therefore isomorphic to the special linear group SL(2,C). In fact, the Mobius 
group of automorphisms of the Riemann sphere is isomorphic to SL(2, C)/Z 2 . There is also an 
infinite set of infinitesimal transformations given by 

£ n : z ^ z + ez n+ \ $ n :z^z + ez n+1 , n G Z. (26) 
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The generators l n = —z n+1 (d/dz) and £ n = —z n+1 (d/dz) satisfy the commutation relations 

[i m , in] = (m - n)i m+n and [l m ,I n ] = (to - n)£ m+n . (27) 

These relations together with the commutator [£ m ,£ n ] = form an infinite-dimensional Lie 
algebra called the Witt algebra or the classical Virasoro algebra. Upon quantization, the gener- 
ators £ n , in become operators L n ,L n , and the commutation relations (27) suffer from ordering 
ambiguities. They require an extension by normal ordering, which occurs for the operators Lq 
and L (when m = —n). The extended algebra is given by the commutation relations 

l)S m ,-n (28) 

l)S m ,-n (29) 

for the generators L n . These extended commutation relations together with the [L m ,L n ] = 
form the quantum Virasoro algebra. The constant c is called the central charge, or conformal 
anomaly. It commutes with all the generators L n and L n . 

Like any quantum field theory, a two-dimensional conformal field theory is determined by the 
space of states ip n which form a Hilbert space , and the collection of the correlation functions 
that are defined within some dense subspace fcW,. The highest weight state of any conformal 
field theory is the state ip whose Lq, Lq eigenvalues are smallest; their eigenvalues L$ip = hip and 
L~otp = hip are the conformal weights, and determine the conformal transformation properties 
of ip. The sum A = h + h is called the scaling dimension, and the difference s = h — h 
is called the planar spin. Let V(ip;z,z) be a field associated to the state ip that satisfies 
ip = y(^;0,0)|0), and |0) is the (SL(2, C)/Z 2 -invariant) "vacuum" state. Then under the 
conformal transformations z — > f(z) and z — » f(z) the field V(ip; z, z) transforms as 

- h 

V(iP;f(z)J(z)). (30) 

A field that transforms in this way is called a primary field. It turns out that there is a 
subspace of states in 7i* that transform as the vacuum state with respect to L n (L n ). These 
are known as chiral (anti-chiral) states. The chiral (anti-chiral) fields then only depend on z 
(z). Therefore chiral fields are given by V(ip; z, z) = V(ip; z) and anti-chiral fields are given by 
V(ip;z,z) = V(ip;z). Correlation functions (n-point functions) of the conformal field theory of 
n primary fields V(ip n ] z n ,z n ) with conformal dimensions h n and h n are represented by 

(0|V(V>i; zi, zi) • • • V{ip n ; z n , z„)|0) = zi, zi) • • • V(ip n ; z n , z„)>. (31) 

The one-point function is given by 

(V{t,z,z))=C (32) 
when h = h with C a constant, and is zero otherwise; the two-point function is given by 

(V(ip i; z 1 ,z 1 )V(iP 2 ;z 2 , z 2 )) = C( Zl - z 2 )- 2h {z 1 - z 2 y 2 ' 1 (33) 



[L m , L n ] = (m - n)L m+n + —m(m 2 - 

for the generators L n , and 

— — — c 

[L m , L n ] = (m - n)L m+ n + T^w(m 2 - 



V(iP;z,z) 



df_ 






dz 




dz 
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when hi = hi and hi = hi, and is zero otherwise. Similar relations can be found for the 
higher n-point functions. In general the Mobius symmetry is used to restrict these functions. A 
remarkable feature of the WZW model is that the correlation functions of chiral fields factorizc 
holomorphically [28]. Let (£ e {1, 2, . . . , N}) be the primary fields of the WZW model 

with conformal weights (h^,h^). If the number N is finite then the conformal field theory is 
said to be rational. In this case the fields are tensor products of two finite-dimensional repre- 
sentations of a semi-simple Lie group with left-hand and right-hand representations labelled by 
and q^, respectively. An ./V-point correlation function of the primary fields is then 

t (V; zi,zi) ■ ■ ■ *%% (V; z N , z N )) = J2 C ab rr- qN ■ ■ ■ > zs)??-™ (z u ..., z N ) 

a,b 

(34) 

with constants C ab that are related to the structure constants of the operator algebra. The func- 
tions T{zi, . . . ,Zn) are the conformal blocks and are completely determined by the conformal 
symmetry. 

2.4 Canonical Quantization of Chern-Simons Theory 

Now, we have seen how the Chern-Simons gauge theory induces the WZW theory on dM.3. The 
relationship is far more profound when this theory is quantized. We now proceed to demonstrate 
this relationship. To this end, let us fix the topology of the three-manifold to A4$ = a x R 1 , 
and let Ao = 0. Here Aq is the component of the connection A in the -Redirection. In this 
gauge the Chern-Simons action (15) becomes 

s ^ = iJ dt J/^{ A 4t A ^ ^ 

and there is one primary constraint 5Scs/5Aq — e^Fij = 0. The next step would be to do a 
Legendre transform to obtain a function that is first order in time derivatives. However, the 
Lagrangian above is already first order in time derivatives and the connections Ai and Aj are 
thus canonically conjugate to each other. Now recall from the definition of the Poisson bracket 
(see Appendix B) that the canonical coordinates {q a ,p a }a=i f° r a discrete classical system with 
a finite number of degrees of freedom will satisfy {q a ,p a } = 1- For a continuous system with 
phase space defined by fields {</>", 7r a }™ =1 , this generalizes to {4> a (x), n a (y)} = S d (x — y) (in d 
dimensions). For the gauge-fixed action (35) we can read off the (equal-time) Poisson brackets 
as 

4-7T 

{A a i {x),A b j {y)} = -e ij 5 ab 5\x-y). (36) 

Now we impose the constraints Fij = 0. This leads to a finite-dimensional phase space given by 
the configuration space of gauge equivalent classes of flat connections. In other words, if T is the 
infinite-dimensional phase space oiA/Q, then the constraint surface F C F is the configuration 
space of Ao/Q- Witten calls this configuration space "a flat bundle on moduli space". We can 
assume a complex structure on a which then amounts to the separation of the connections Ai 
into holomorphic and anti-holomorphic degrees of freedom, as was done above for (18). The 
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configuration space Aq /Q is then a flat vector bundle on moduli space. It turns out that these 
bundles are exactly those that arise naturally in the holomorphic factorization of the Appoint 
functions that we have already defined above as the conformal blocks T(z\, . . . , zn)- We now 
see the deeper relationship between a generally covariant gauge theory in three dimensions and 
a conformally invariant quantum field theory in two dimensions: quantization of the Chern- 
Simons action on M3 = crxR 1 produces a Hilbert space Ho of wavefunctions that are recognized 
as the generating functional for two-dimensional current correlator blocks. 

3 The Theory of Knots 
3.1 Topological Preliminaries 

A knot is a submanifold of K 3 that is diffcomorphic to the circle S 1 . The simplest example of 
a knot is the unknot, a circle that contains zero crossings: 

{(x,y,z)GM. 3 \x 2 +y 2 = 1,^ = 0} ~ O- ( 3? ) 

Another example of a knot is the trefoil knot ■ This is an example of an alternating knot 

where the crossings alternate between over and under. A link is a submanifold of M 3 that is 
diffcomorphic to a disjoint union of circles. The circles are the components of the link. A knot 
is a link with one component. An example of a link is the Hopf link (JT) . Two knots K 
and K' are said to be ambient isotopic if! K can be continuously deformed into K' through 
embeddings of R 3 . An equivalence class of links with ambient isotopy as the equivalence relation 
is said to be an isotopy class. A knot is amphicheiral iff it is isotopic to its mirror image. The 
trefoil, for example, is not amphicheiral. Also, a knot decomposes if it is the connected sum of 
two knots. A knot K is said to be the connected sum of knots K\ and K 2 , denoted Ki#K 2 , if 
a two-sphere can be embedded in M 3 so that the two-sphere intersects K in exactly two points 
so that K is bisected into pieces one of which is K\ minus a small arc and K 2 minus a small 
arc. A knot is then prime if it cannot be decomposed into two simpler knots. 

It turns out that two knots K\ and K 2 are isotopic iff there exists a sequence of elemetary 
isotopies (called Reidemeister moves) between projection diagrams of K\ and K 2 . The Reide- 
meister moves are given in diagrammatic representation in Fig. 1. The first move RI consists 
of modifying a diagram by putting a twist in the neighbourhood of a single strand. The second 
move RII consists of cancelling a pair of crossings. The third move consists of sliding a strand 
over a crossing (RHIa) or under a crossing (RHIb). There is also a "zeroth" move R0 which 
consists simply of an isotopy of the plane. We note that these moves are local; they are applied 
(one at a time) to one segment of a knot diagram while leaving the rest of the knot diagram 
unchanged. The equivalence relation generated by all Reidemeister moves except RI is called 
regular isotopy. Also, a loop that cannot be deformed into the unknot with the Reidemeister 
moves is said to be knotted. 
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Figure 1: The Reidcmeister moves RI-RIII(a,b). Here a strand is a single line segment. 
3.2 Link Invariants 

By Reidemeister's theorem [29], the topological properties of knots and links in three-space (M 3 
or S 3 ) can be studied by working with diagrams of the knots taken up to the Reidemeister 
moves. That is, if diagrams for two knots or links are ambient isotopic in three-dimensional 
space iff the diagrams can be obtained one from another by a sequence of Reidcmeister moves. 

A knot or link is said to be oriented if a direction along the curve is chosen for each component 
of the knot or link. When we have a diagram of an oriented knot or link, then each crossing 
inherits and orientation sign. We define for right-handed orientation positive sign e( ^ ) = +1 

and for left-handed orientation negative sign e( ^ ) = — 1. Let us first consider the self- linking 
number. For an oriented (one-component) knot diagram K with TV crossings {1,2,..., N} this 
is 

N 

w(K) = ]Te 9 . (38) 

9=1 

We note that this function is not preserved under RI but is preserved by the other Reidemeister 
moves, and is therefore a regular isotopy invariant. Now let us consider an oriented diagram 
formed by two knots K\ and K 2 . The linking number is 

L(K U K 2 ) = \ ( 39 ) 
p£K 1 nK 2 

where the intersection KiC\K 2 does not include the self-crossings of K\ or K 2 . As an example 
we will calculate the linking number for Hopf links with opposite orientations. For K\ ~ 
and K 2 ~ so that K\ U K 2 <~ (j±T) we have 

L(K u K 2 )= 1 -(e(X ) + e(X ))=+!, (40) 
while for K x ~ >Q and K 2 ~ {J so that K 1 l>K 2 ~ (£f) we have 

L{K u K 2 ) = \(eCA ) + <"A ))=-!■ ( 41 ) 

Before we go on to define the bracket polynomial we will introduce the notion of framing. 
Let if be a knot that is defined by a closed curve C. The framing of K is a vector field that 
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Figure 2: The modified reidemeister move RI' for a framed knot segment. 

is normal to the curve C, whose endpoints generate a boundary curve C . So a framed knot 
can be visualized as a closed ribbon with boundaries defined by the curves C and C . The first 
Reidemeister move RI fails for framed knots. Instead, a curl in a diagram corresponds to a 2tt 
twist in the ribbon. See Fig. 2. We find that both [a] and [b] represent the same 2ir twist and 
so we set [a] ~ [b]. This modified deformation is henceforth labelled RI'. 

Now we are ready to define the bracket polynomial. Let £ be a multi-component link 
containing N crossings. A state of C is obtained by choosing a smoothing of each crossing. 
Assigned to each crossing of £ is a choice a q £ {A, B}, where the action of the vertex weights 
A, B on the gth crossing is given by 

X - A(X) and X - £()(); ( 42 ) 

smoothing the crossing horizontally results in a multiplicative factor of A while smoothing the 
crossing vertically results in a multiplicative factor of B. The resulting state will be a finite set 
of disjoint circles of cardinality a. The bracket polynomial is then a state sum given by 

( £ )=E d " i n^- ( 43 ) 

a q 

This is a polynomial in the commuting variables A, B and d. The skein relations that define 
the bracket polynomial are 

(X)=MX)+ B and <O u K)=d{K). (44) 

These equations follow directly from the definition of the state sum, where a normalization of 
unity has been employed for the unknot. Here K does not intersect the added loop. Note that 
the bracket polynomial is an invariant of regular isotopy; if the axioms are applied to a twist 
then 

(D-) = (Bd + A){ — ) and (£>_) = (Ad + B)( — ) . (45) 

In fact, the bracket polynomial is invariant under the framed Reidemeister moves for a specific 
choice of the commuting variables A, B, d. As defined, the bracket polynomial is invariant under 
RI', but fails for RII: 

(JT)=AB() ()+[A 2 + ABd + B 2 ](X)- (46) 
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Invariance requires that AB = 1 and A 2 + ABd + B 2 = which uniquely determines the 
variables such that B = A~ x and d = — (A 2 + A~ 2 ). With this choice the bracket polynomial 
is now also invariant under the action of RIII. Now we can write the defining state sum (43) 
for framed links that is invariant under the Reidemeister moves as a polynomial in the single 
variable A, with the skein relations 

(X)= A (X)+ A ~ 1 and (O U K)=-(A 2 +A- 2 )(K). (47) 

Let us now define the Jones polynomial. For a framed oriented link C this is a polynomial in 
A such that 

V C (A) = (-A- 3 ) w( - c \ C )(A). (48) 
The skein relations for this polynomial (with q = A~ 4 ) are 

?X -T'X =(^-T = )X and O =1 - (49) 

Now note the relation for the bracket polynomial under the action of RI for unframed links 
given by (45), with coefficients Bd + A = -A~ 3 and Ad + B = -A 3 . Under RI this extra 
factor acquired by the bracket polynomial is cancelled by the coefficient in the definition of the 
Jones polynomial above. The Jones polynomial is therefore an invariant of oriented links. The 
framing is irrelevant. 



3.3 Link Invariants from Gauge Theory 

The Jones polynomial was originally defined for knots and links in the three-sphere S 3 . We 
will proceed to show that the exponential of the Chern-Simons action as a measure provides 
an intrinsic definition of the Jones polynomial in an arbitrary three-dimensional manifold Ai^. 
The reader should note that the functional integral so defined is only understood heuristically. 
It is nevertheless, a powerful heuristic device, leading to a host of results. We might expect 
this because the Chern-Simons action is invariant under diffcomorphisms, so a change of basis 
to loop space (see below, where the loop transform is defined) should give a knot invariant for 
every gauge group G. Let C be an oriented closed curve in M3. We consider the Wilson loop 



Wr[C] = Tr fl Pexp 



Aidx 1 

c 



(50) 



obtained by calculating the holonomy of the connection Ai in the representation R of the gauge 
group G. Clearly the closed curve C can be a knot; taking r oriented and non-intersecting knots 
C| (£ € {1, 2, . . . , r}) whose union is a link £, we define the partition function of the manifold 
M3 such that 

Z{C)= [ (DA)exp(iScs)f[WR i (C 6 ). (51) 
JA/g 



This is a path integral over all gauge orbits where we have assigned a representation of G to 
each knot Q . The partition function above is the unnormalized expectation value of the given 
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Figure 3: Cutting of the manifold S 3 with Wilson loops into two manifolds Ml and Mr. 



Wilson loop. Evaluating this path integral amounts to the determination of a vector tp in the 
physical Hilbcrt space Hm 3 , i-e. the space of gauge-invariant wavefunctions in the presence of 
Wilson loops. We will now consider the special case My, = S 3 and G = SU{N). Assume that 
there is a link present in S 3 and we cut S 3 into two sections Ml and Mr. The boundaries 
of these sections will now be Riemann surfaces intersecting the link transversely in a set of 
marked points. In the case where there are four marked points and the Wilson lines are all 
in the same TV-dimensional representation of SU(N), then (we state without proof that) the 
physical Hilbcrt space will be two-dimensional. If the sections Ml and Mr are pasted back 
together, then the resulting link invariant may be the original, or may be a result of pasting 
Ml with X\ or with A2. This is illustrated in Fig. 3. 

Now let ipL € Hm l and ipR, ipi,ip2 & 1~Lm r - Here tpR is a vector associated to the boundary 
of Mr, ip\ is a vector associated to the boundary of X\ and V2 is a vector associated to the 
boundary of X 2 . We note that any three vectors in a two-dimensional vector space are linearly 
dependent, and since ipR,ipi,ip2 & T^Mr are vectors contained in the same Hilbert space we 
have that 

aipR + (3i>! + 7^2 = 0. (52) 

Here the constants a, /3, 7 G C . We may then postulate a natural pairing of vectors (ipL,ipR) 
to the partition function Z(C), (ipL,4>i) to the partition function Z{C\) and (if>L,4>2) to the 
partition function Z(£ 2 )- The relation (52) becomes 

olWlM + 0W>L,ih.) + 7(^,^2) - 0, (53) 

and this translates into a reccursion relation for the expectation values of the Wilson loops 
given by 

aZ{C) + f3Z{C 1 )+- / Z{C 2 ) = 0. (54) 
Diagrammatically this is interpreted as 

«X+^+7%=0- (55) 

Now, the explicit evaluation of the constants a, f3, 7 would divert us too much from the main 
purpose of this paper, but it is neccessary to state the final result. The reader is referred 
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to Witten's paper [3] and references therein for details of the evaluation. For a connection 
A of Chern-Simons theory with gauge group G = SU(N) and all Wilson lines in the same 
TV-dimensional representation of SU(N), we have that 

( 2?ri \ 

a = — exp — — — — I (56) 



\N{N + k) 
fnr(2-N-N 2 )\ ( m{2 + N - N 2 ) \ 

= ~ eX H N(N + k) j +eX H N(N + k) ) (57) 

7 = eXP V N(N + k) J' (58) 

Multiply these by a factor exp[jri(N 2 -2)/N(N + k)] and define a variable q = cxp[2iri / (N + k)} . 
Now (54) and (55) become 

-q^Z{C) + {q^-q-^)Z{Cx)+q-^Z{C2) = (59) 

This is the skein relation for the HOMFLY polynomial [30], which is the generalization of 
the Jones polynomial to G = SU(N). If G = SU{2) we get the skein relation for the Jones 
polynomial. 

Before moving on to gravitation, we will briefly discuss a simple case of the partition function 
(51) in order to illustrate a subtlety that arises when evaluating the path integral (51) for a 
specific case. Take G = U(l) and Ms = S 3 . The potentials commute so the cubic term AaAaA 
in the Chern-Simons action vanishes, while the trace and path-ordering become irrelevant. 
Consider in 5* 3 two non-intersecting knots C a and Cb- To the two different representations of 
U(l) we assign integers n a and rib- Evaluating the integral along the knots C a and Cb in a 
region of S 3 that is essentially Euclidean space with x % and yi the coordinates along the knots, 
the partition function for the case a ^ b becomes 

Z(C a ,C b ) = ±-[ dx l [ dy^ t3k ^—^. (61) 
4fc7c a Jc b \x~y\ 3 

This is the Gauss linking number, which was first written down by Gauss more than one hundred 
and fifty years ago! It is well defined provided that the two knots do not touch or intersect, 
because the integral diverges at the point x = y. Witten regularized this divergence for the 
case a = b by giving the knots a framing, and interpreted the resulting quantity as the self- 
linking number of one knot. Actually, this framing regularization is neccessary for all gauge 
groups. We have not mentioned it before because it was not relevant to the derivation of the 
skein relation, i.e. we did not need to evaluate any integrals. In a sense, we found a rcccursion 
relation for Wilson loop expectations for unframed loops, thereby neglecting the self- linking. 
That is why we obtained the Jones Polynomial. Later we will see that the need for framing will 
arise naturally for a state of quantum gravity with non-zero cosmological constant. There we 
will obtain the bracket polynomial for G = SU(2), which depends on the framing. 
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4 The Gravitational Field: Classical Theory 



4.1 Second Order Formulation 

The use of tensors in general relativity is motivated by the principle of general covariance 
[31,32], which states that the laws of physics must be independent of the coordinates being 
used. For example, the metric tensor relates distances in space and time via the invariant 
line interval ds 2 = g )lv dx il dx v ', and may be seen as a generalization of the Pythagorean theorem 
to curved four-dimensional spacetimc. The metric tensor is normally taken to be the dynamical 
field of the gravitational force. The equations relating the gravitational field to its source are 
the Einstein equations 

R^u - \R9» V - Ag^ = 8irG N T^. (62) 

Here the Ricci tensor R^ v encodes information about the curvature, and its contraction R = 
g^Rfiv is the Ricci scalar. The material contents of the system, whatever it might be, are 
described by the stress-energy tensor T^. The quantity A is the cosmological constant, which 
manifests itself as an acceleration whose magnitude increases with distance r as |A|r/3 in the 
weak-field limit. The dimensions carried by A are L~ 2 , which implies that this term must 
be significant globally. The Einstein equations (62) can be derived from an action, either the 
Einstcin-Hilbcrt action 

$ matter 

(63) 

by taking the Ricci scalar to be a functional of the metric, or equivalently the Hilbcrt-Palatini 
action 



5[5 ^' r ^ J = 16^7 / d V |det(5 ^ )l(i?[r V] - 2A ) + S ' 



matter 



(64) 



by taking the Ricci scalar to be a functional of the Christoffel coefficients of the second kind 
r 7 ^^. The matter action S ma tter — SmatteAg^,^^,^] is generally a function of the scalar 
field (j), spinor field tp and vector field A in addition to the metric g^ v . The stress-energy tensor 
is obtained by varying the matter action with respect to the metric such that 



2 5S n 
^~5g^ 



T, v = ^T,T r - (65) 



It must be noted, however, that this tensor is not the same as the canonical stress-energy 
tensor that is obtained from the symmetries of the Lagrangian, i.e. Noether's theorem. The 
equivalence only holds for a scalar field (up to a numerical factor). In general, the canonical 
stress-energy tensor will not even be symmetric, and in many cases such as electrodynamics, will 
not even be gauge invariant. Thus we define the stress-energy tensor as that which is obtained 
by variation of the matter action with respect to the gravitational field. When the action (63) 
and the action (64) are varied with respect to g^ we get the Einstein equations. The virtue 
of using (64) over (63), however, is that variation of the former with respect to V 1 also gives 
the metric compatability condition V 7 <7 M „ = 0. 
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4.2 First Order Formulation 



The Einstein theory can also be formulated in terms of an orthonormal tetrad frame ej that 
satisfies 

S^e/e/ = f] IJ and ??/je/e/ = g^. (66) 

Locally then the manifold is flat Minkowski space M 4 and the curvature is encoded in the tetrad 
frame field. The covariant derivative acts on spacetime indices in the usual way: 

= d^v x + T\ v v l/ . (67) 
In addition now the covariant derivative acts on Lorentz indices such that 

= d^v' +uj/jv j . (68) 

This defines the spin connection w^ 1 j. The spin connection determines the torsion two- form 

T 1 = de 1 + lo'j A e J = 0; (69) 

the vanishing of the torsion implies the antisymmetry oj ij — —oj ji of the spin connection. This 
with the orthonormality of the tetrad eje„i = i]^ leads to the metric compatibility condition 
V 7 <7 M1 , = and determines a unique connection on the manifold called the affine connection. 
With this connection are associated the connection coefficients r . In terms of these the spin 
connection coefficients are 7 = Y 1 JK e^ K . This is general and may be taken as the definition 
of the r-coefficients. In a coordinate basis these are the Christoffel symbols of the second kind. 
Now it is clear that the metric compatibility condition implies that the spin connection is also 
compatible with the metric: 

V„e/ = d v ej - V^eJ + Lojje/ = 0. (70) 

This can be inverted to give back (69). The spin connection also determines a curvature 

V J = W J - 9»«>" + """J J - (71) 
the components of which determine the curvature two-form 

R 1 j = Rj jdx^ A dx v = du 1 ' j + oj'k A lu K j. (72) 

Equations (69) and (72) are called the Cartan structure equations. If the tetrad is invertiblc 
then we can define the Ricci tensor Rj = R^J" 1 e" ' j and Ricci scalar R = Rje^j. With this 
the Hilbert-Palatini action becomes 

S[e,u;} = —±—[ d i x\det(e)\{e^ I e\,R^ IJ - 2A) + S matter , (73) 

167TGAT J M 

which up to boundary terms is invariant under local SO (3, 1) transformations and local trans- 
lations. The matter action S matter = S matter [e,<j>,ip,A] is generally a function of the scalar 
field (j>, spinor field and vector field A in addition to the gravitational field e. Variation of 
the total action with respect to e 1 gives the Einstein equations 

e"j V" - l^Me P N e v J R w MN - e/A = 8nG N T„ 3 (74) 

as the equations of motion, while variation with respect to uj ij gives the metric compatability 
condition. Here the stress-energy tensor T^ 1 = SSmatter/Se^j defines the three-form T 1 = 
T fl I e tl vpX dx v AdxP Adx x . 
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4.3 Self-Dual Formulation 



We finally come to the self-dual formulation of general relativity [33-35] , which is the starting 
point for canonical quantization. To this end, we consider the integral 

T[eM= I d i x\dct(e)\e^ I er j e IJ MN R^ MN (75) 

which is a topological term in the sense that it is invariant under local variations of the tetrad. 
An equivalent action for general relativity is therefore 



+ S[eM = S [e,u:}-^-T[eM 



1 



d 4 a;|det(e)|eV / 



D IJ 1 J J D MN 



IQttGn Jm 

Now in terms of the spin connection uj the self-dual connection is given by 

1 



1 - 2 



I J " I J MN 



J H 2 MN^fi 

and this determines a self-dual curvature 



(76) 



(77) 



F IJ - r) A I J - r) A I J 4- A IM A J - A IM A J 

— E> IJ 1 JJ TJ MN ( 7 o\ 

With this and a change of variable uj — > A the self-dual action becomes 

+S[e,A] = J d*x\det(e)\(e' 1 I e''jF llv IJ - A) + S matter . (79) 

Variation of this action with respect to the tetrad field does give the Einstein equations as can 
be verified. Here the connection is complex so the metrics that satisfy the Einstein equations are 
complex- valued. To recover real general relativity certain reality conditions must be imposed 
on the phase space. In particular, we begin with complex general relativity on a real manifold, 
and, after defining the Hamiltonian using the self-dual connection, we restrict ourselves to the 
real section of the complex phase space. This is not very straightforward in practice, however, 
because it is not obvious that the pull-back of the symplectic structure to the real section is 
itself real. An alternative approach is to start with a real Ashtckar connection for Euclidean 
general relativity and later Wick rotate to Lorentzian general relativity. A more modern ap- 
proach, however, is to use a form of Ashtekar variables with a real connection that describes 
Lorentzian general relativity, but at the price of having to solve a much more complicated 
Hamiltonian constraint. We will proceed to derive the gravitational Hamiltonian using this 
Barbero formulation [36] that describes the phase space of Lorentzian general relativity with a 
real connection. Later in Section 5.3 we will go back to the complex self-dual connection and 
derive the Kodama state for the simplified Hamiltonian constraint. 



4.4 Classical Phase Space 

The canonical structure is obtained using the usual ADM decomposition (see Appendix C). 
Without loss of generality we use the internal gauge symmetry of the gravitational field to fix 
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the tetrad component e° - = 0. Now define on the time-slice £ the densitized triad E l - — ^fqe % - 
of weight one and local SO(3) connection = e 0IJK A i j^, which are related to the phase 
space of the induced metric and extrinsic curvature by 



PjEP^qq* and Aj = T i i + iK^ . 



(80) 



Here T i I is the SO(3) connection of the triad. From this we can consider a more general 50(3) 
connection given by [36] 



V = r l i + 1 K l i , 



(81) 



where 7 is an arbitrary parameter. If 7 is real it is known as the Immirzi parameter. It does 
not play any role in classical dynamics because the extrinsic curvature term does not show up 
in the equations of motion unless the action contains fermions. (It was shown in [37] that 7 
appears in the equations of motion as the coupling constant of a four-point fcrmion interaction 
term.) In terms of the densitized triad E % . and the connection in (81) the action + S becomes 



s = 



1 



dt / d 6 x 



1 N 



Here we have seven first class constraints 



1A/E 1 f - iA oi g' + iN% - -—S 
1 2^ 



Q 1 = 
V, = 



D t E u ft 
E 3 .F- 1 

I l 3 











S = e Uk E i i E i 



2* 



F- 



ijK 



-^SijkF k 



2(1 +7 2 )^ 



& [i E^ ] (A l I -T l I )(A/-T J J )^0, 



(82) 



(83) 
(84) 



(85) 



where Di and are the gauge covariant derivative and curvature of the connection A i 1 . (The 
notation "w" denotes weak equality on the constraint surface.) In the above action we identify 
the term iA i 1 E % - as the analogue of the term in the Legendre transform / "p^cf — H. 

Thus the gravitational (Ashtekar) Hamiltonian is a linear combination of the constraints 



dt / (Fx 



1 TV 



-.S 



(86) 



The phase space is now described naturally by the SO(3) connection A i I (gauge potential) and 
its conjugate momentum E % ~ (electric field). They satisfy the Poisson bracket relations 



{Aj,E j j} - ^GnS^IS^x-x') 



(87) 
(88) 



Now we know that the first class constraints of a Hamiltonian system generate gauge trans- 
formations. General relativity is no exception. The three vector (diffeomorphism) constraints 
Vi and scalar (Hamiltonian) constraint S generate surface deformations that are equivalent to 
spacetime diffeomorphisms when the field equations hold. The algebra of these surface defor- 
mations is a Poisson algebra with structure functions rather than structure constants. We have 
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now, in addition, the three "Gauss law" constraints Q 1 that are typical of Yang-Mills theory. 
These generate the usual 5*0(3) gauge transformations. 

We make a final note on the canonical theory described here. The formulation in terms of 
the connection A i 1 turns out to be a generalization of the Einstein theory. The equations in the 
metric formulation contain the Ricci scalar which requires the inverse metric to be defined. This 
means that the metric has to be non-degenerate. The phase space here is now determined by 
the set (A i I , E l .), so the metric is a derived quantity. Here all the equations of the Hamiltonian 
theory are polynomial, and the inverse triad E { 1 does not appear in any of them. The theory 
is therefore well-defined even when the triad is degenerate, and therefore the theory admits 
solutions with a degenerate metric. If the triad is restricted to be non-degererate then the 
Einstein theory is recovered. 



5 The Gravitational Field: Quantum Theory 
5.1 Loop Representation 

The phase space of general relativity is now defined and the canonical analysis is complete. 
What is left is to quantize the system. This is done by promoting the constraints (85) to 
operators Q 1 — > Q 1 , Vi — > V«, S — > S: these are required to annihilate physical states on the 
Hilbcrt space. In the connection representation the phase space variables become 

A/^A- 1 and & f -» & t = -8njG N (89) 

Also the Poisson brackets on the phase space are promoted to commutators on the Hilbcrt 
space via {x,y} — ► (l/ih)[x,y\. Of course this is quite difficult in practice due to ordering 
ambiguities. General relativity also suffers from the problem of time, structure functions of the 
surface deformation algebra etc. We will not go into the details here, but refer the reader to 
the book by Ashtekar [35], the review article by Rovelli [38] and the book by Rovelli [39]. 

The quantum constraints require that states be gauge invariant functionals of the connection. 
And so we are naturally led to consider Wilson loops defined by the self-dual connection. 
To this end we consider an SU(2) spinor description for the canonical variables where Ai = 
A i jT I and E % — 2E % -t 1 in terms of the spin half SU(2) generators t 1 = —(i/2)<r I in the 
fundamental representation. (Here a 1 are the usual Pauli matrices.) We consider along a 
continuous piecewise smooth curve 7 : [0, 1] — > £ that is parametrized by s — ► {7 a (s)} (a e 
{1,2,3}), the holonomy of the Ashtekar connection Ai such that 



U 1 (s 1 ,s 2 ) = Pexp 



This satisfies the differential equation 



S2 



. dx l (s) . j 
ds 



(90) 



j-U 7 (s, s ) + ^Ma^P^s, s ) = (91) 

and the boundary condition J7 7 (s ,so)a = ^f; this simply says that the vector U 7 (s) in the 
vector bundle over the manifold M is parallel transported along the curve 7. Hence the alternate 
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name "parallel propagator" for the holonomy. Now we define the fundamental set of gauge- 
invariant loop operators (for A = 0) [38,39]. For a loop 7 such that 7 l (0) = 7*(1) and the 
points si, s 2 , ■ ■ ■ , sat <G 7 we have: 



T[ 7 ] = -Tr[*7 7 (0,l)] 
T%]( S ) = -It [l7>, 



T^K*!,^) 



-Tr 



f/ 7 (s 1 ,s 2 )^( S2 )f/ 7 (s 2 ,s 1 )^ 1 (si) 



T ll - l «[ 7 ]( Sl ,..., SJ v) = -Tr U 1 (s 1 ,s N )E^(s N )...U 7 (s 2 ,s 1 )E i ^s 1 ) 



(92) 
(93) 

(94) 

(95) 



We will now proceed by introducing some notation for the loop variables that will be needed 
below. Let a and (3 be two loops, and I be a segment of a curve. If a and (3 intersect at a point 
p then the loop ajj/3 is obtained when starting at p, going first through a, then through /3 and 
finally ending at p. If a segment of a curve £ is glued to the loop a then the result is denoted 
a o I. Also the inverse of a loop or 1 or of a segment of a curve i" 1 is just the loop or segment 
with opposite orientation. With this, the loop operators satisfy the properties: 



T[a] = Tla- 1 } 

T[a] = TfaoW- 1 ] 

T[a]T[(3] = rwi + rnr 1 ] 

limTfaf.1 = 2. 



(96) 
(97) 
(98) 
(99) 



Here a% is a loop centered at the point x on the spatial manifold £ with area e. Finally we 
write down the (closed) Poisson algebra of the lowest-order operators T and T % : 



{T[a},T[/3}} 
{T>]( s ),T[/3]} 







{T>](s),T>](i)} = --A%a{s)] [Ti{a$ s m)-'r i ns[3- 1 }(t)} 



(100) 
(101) 

(102) 

(103) 



Here we have introduced the notation jj s to mean that the breaking and rejoining of the loops 
occurs at the intersection of the two loops where the parameter is s, and we have defined the 
structure functions 



A l [f3,x} = j> dt/3 l (t)6 3 ((3(t),x). 



(104) 



Note that if the loops do not intersect then the Poisson brackets vanish. These can be obtained 
with rigor by using the operator definition for the E l ~ field. It is to be understood that these 
relations define the Poisson brackets and are not the usual Poisson brackets that are defined by 
derivations. The relations above are called the T-algebra, and realizes Isham's proposal [13] of 
employing a non-canonical Poisson algebra for the quantization of general relativity 

Let us now prove that loop states solve the constraints of the phase space of general relativity. 
We will do this for any arbitrary loop state. We have already encountered a gauge invariant loop, 
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namely the Wilson loop, which is just the trace of the holonomy of the Ashtekar connection. 
The gauge invariance implies that this quantity automatically solves the Gauss constraint Q 1 . 
In addition, it is an exact solution to the scalar constraint S. So let us consider a general 
functional of the form 



H[{~/},A] =H TrPexp I A 



7£{ 7 } 



(105) 



where 7 is a loop contained in the set of differentiable non-intersecting loops {7}. Then for a 
functional &[A] in the connection representation and a functional ^[{7}] in the loop represen- 
tation, we define the loop transform [10, 11] 

*[{7}]=/dM^[{7},A]*[A]. (106) 

Here d/z[A] is a measure on the space of connections that is assumed to be invariant under 
diffeomorphisms. The transform is therefore a mapping T : $ — > ^[{7}] = &[H [{7}, A]]. 
Operators in the two representations are related to each other through the transform T in the 
following way. Let O be an operator acting on the space of connections and let O be an operator 
acting on the space of loops. Then O and O are related through the transformation 

OT = TO^ (107) 

provided that OH [{7}, A] = OH[{j}, A}. Now, the loop transform is the key to understanding 
how the loop functionals solve all constraints in the loop representation. It is known that the 
Wilson loop satisfies Q 1 and S. Therefore the kernel H[{-f}, A] also satisfies Q 1 and S, and it 
follows that the functional ^[{7}] in the loop representation satisfies the constraints for any 
choice of functional &[A] in the connection representation. All that remains now is to impose 
the vector constraints V». However, diffeomorphisms map loops to other loops in the same 
knot class, so we conclude that ^[{7}] will satisfy all constraints provided that the functional 
depends on the knot classes of loops. 



5.2 The Spin Network Basis 

The most beautiful result of loop quantum gravity is the discretization of the spatial manifold. 
The non-perturbative quantization of the gravitational field leads to quantum geometry oper- 
ators with discrete spectra. Details on this can be found in [39,40] and references therein. A 
self-contained exposition on the calculation of geometry eigenvalues based on Temperley-Lieb 
recoupling theory [41] can be found in [42]. The area operator, for example, has been used to 
calculate the Bckcnstcin-Hawking entropy for an extremal black hole [43-46]. This result was 
used to fix the Immirzi parameter. 

Functionals in the loop representation that diagonalize the area and volume operators are 
represented by combinatorial structures called spin networks. A spin network state [42] is a 
combination of loops called a graph T ~ {7}, containing a certain number of edges and vertices. 
To each edge we assign a half-integer (spin) ji which labels an irreducible representation of 
SU(2), and to each vertex we assign an invariant tensor v a called an intertwiner in the tensor 
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product of the representations ji, . ■ . , j n at which the spins join. A spin network is then the set 

S = {T,j it v a }. (108) 

Note that S is an abstract spin network. The embedding of S into the spatial manifold is 
an embedded spin network. The distinction should be clear from context. Essentially we 
find the holonomy of the connection A in the representation ji and use the tensor v a at each 
vertex to obtain an invariant functional. The result is then a state ^[A] in the connection 
representation. The spin network states are linearly independent and provide an orthogonal 
basis for the kinematical Hilbert space [38, 39]. In fact, we can write any state in the connection 
representation as a superposition \ip) = ^2 s (S\ip) \S). The coefficient (S\ip) here is the inner 
product 

(Sty) = [_d^[AWsW[Al (109) 

J A/Q 

where the diffeomorphism invariant measure d/x[A] is over the gauge equivalent classes of gener- 
alized connections A. The configuration space A/Q is a canonical enlargement of A/Q [47-50]. 
The enlargement is neccessary if the loop operators are to remain well defined in the quantum 
theory where the number of degrees of freedom is infinite. This measure has only been defined 
rigorously for A = 0. In the next section we will discuss the case A ^ 0. There, we will see 
how the loop transform is related to Witten's partition function when we discuss the Kodama 
wavefunction, and obtain a known topological invariant as an exact quantum state in the loop 
representation. As we will see, the presence of A can be incorporated into the measure of 
Chern-Simons theory. 



5.3 An Exact State 

Recall the scalar constraint S in (85). If we set 7 = i as for the self-dual connection, then the 
last term vanishes and we get 



S = e IJK E l jE j j 



rr -L—*F k 

r ijK "t" 3 fc ijfc & k 



0. (110) 



This is the Hamiltonian constraint in the original self-dual formulation. It is clear that a solution 
to this S (and also to Q 1 and V%) will have to satisfy 



If A { f = iy/K/'ifS^ is chosen for the connection with arbitrary function / = f(t), then the 
curvature becomes F^g = — ^ Z 2 ^ Now the connection defined in (80) is purely imaginary 
so the SO(3) connection r f 1 of the triad vanishes which means that S is a flat manifold. The 
canonical momentum becomes E l ~ — f 2 8 l ~ so that the spatial metric is = f 2 5ij. To find 
the evolution equation for / we write down the Hamilton equations of motion for the canonical 
variables: 

V = {V.#a} and fr t = {& t ,H A }. (112) 
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For simplicity we take the shift vector N 1 in (86) to be zero. In other words, the spacetime 
metric has no t-x % cross terms. The evolution for / is then given by / = N \J 'A/3/ 4 , and taking 
the lapse function TV = l/^/q = l// 3 leads to / = ^/A/3/ or / = exp[^/A/3t]. So the line 
element is 



ds 2 = —dt 2 + q i jdx z dx : ' = —dt 2 + exp 



(113) 



which is immediately recognized as de Sitter spacetime. This is also a Fricdmann-Robertson- 
Walker spacetime with flat spatial manifold and scale factor R(t) — exp[y / A/3t]. It is the 
only self-dual solution to the Einstein equations for Lorentzian manifolds. Now, to see how 
this solution is related to topological field theory in S, we take S(A) to be a Hamilton-Jacobi 
functional on the space of connections A in S. This implies that the canonical momentum is 



E l 



-8ttG 



N " 



6S(A) 



but from the curvature (111) we have then 

_ 8ttAG n 

* ' jkl — 3 ei 3k- 

In terms of A, 1 the curvature is 



SS(A) 

saJ ' 



(114) 
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Fij 1 — ^iAj 1 djAf 1 + e^ 1 A"Aj. 
Integration then gives the Hamilton-Jacobi functional 



From this we arc led to the semiclassical state 



^k(A) = A/exp 



8ttAG 



A A dA + -A A A A A 



Y CS 



N Js 



(116) 



(117) 



which describes de Sitter spacetime. This state is, however, an exact state in quantum gravity. 
It is annihilated by all the constraints: 



g^xiA) = 0, V t ^ K (A) = and S^ K (A) = 0. 



(119) 



This means that ^ k{A) is a physical state of non-perturbative quantum gravity with a non- 
zero cosmological constant [51]. This is the Kodama state, whose classical limit is de Sitter 
spacetime. States of the form 



= Wexp 



8ttAGn Jj2 



Jjy cs + 8ttAG n S'[A}) 



(120) 



can then be used to obtain linearized gravity on de Sitter spacetime, i.e. gravitons propagating 
over a classical background. The normalization factor N satisfies 8N /SA — 0, and depends on 
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the topology of the manifold. In fact, Soo [52] pointed out that a more general solution would 
sum over all inequivalent topologies of £ so that 



* [A]= ^L^ exp l^L Ycs 



{Top(£) |<9E=0} 



(121) 



The boundary of £ is required to vanish so that ^f[A] is still annihilated by all the constraints. 
This sum, however, is very large and in fact unknown, because a complete classification of 
three-manifolds is lacking at this time. 

Finally we come to a very striking property of the Kodama state, which serves as a point 
of convergence for all the material that has been discussed up to this point. So let us take the 
loop transform of the Kodama state. Allowing for either Euclidean or Lorentzian manifold M 
by introducing the integer e = {i, 1}, we have 



*[{7}] = / duMHfalAWexp 



* Yes 



8ttAG 



N 



(122) 



We want to be able to evaluate the loop transform for states with A ^ so as to include the 
Kodama state. The above integral, however, is undefined unless the loops are framed. With 
this framing regularization the group SU (2) of spin networks becomes "quantum deformed" to 
SU q {2) with deformation parameter 



q = exp 



2tt 



k + 2 



with fc =Gp/ (123) 



Clearly then SU q {2) — > SU(2) for A — > 0, q — > 1. For Euclidean manifold e = i this is exactly 
Witten's partition function with gauge group G = SU{2) and the Wilson loops all in the 
same two-dimensional fundamental representation of SU(2). In loop space the Euclideanized 
Kodama state is then the bracket polynomial in the variable A. Unfortunately, this picture 
is a bit more difficult for the Lorentzian case e = 1 because now the integral (122) has to be 
taken over a contour since the connection is complex. It is possible to analytically continue 
the deformation parameter to complex values of the Chern-Simons level. With the loops now 
framed, the spin networks become essentially blown up to a two-dimensional surface [53-55]. 
Edges become tubes with ruling. The ruling means that the original loop before deformation can 
be identified. Vertices become spheres with punctures, each puncture being the location where 
a tube intersects the sphere. The punctures are now labelled by intertwiners of the quantum 
group G q = SU q (2). The intertwiners depend on the level of the Chern-Simons theory, and are 
none other than the space of conformal blocks. The reader will find this easy to visualize if 
he/she thinks of the spin networks as just trivalent graphs and replaces the edges by tubes and 
the vertices by junctures of tubes that form non-singular surfaces. We refer the reader to the 
literature cited in this subsection for the handling of framing with respect to the networks and 
the surfaces. 

Now the picture is complete. Functionals of non-perturbative quantum gravity with a pos- 
itive cosmological constant in the loop representation are finite-genus ruled Riemann surfaces 
with marked points. It is remarkable that the structures which naturally describe perturba- 
tive string theory [19, 56] have emerged as functionals of a non-perturbative quantization of 
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the gravitational field. One thing still lacking in perturbative string theory is a background- 
independent formulation called M(ystery) theory. An ambitious program launched by Smolin 
[57-60] is to describe M theory as a deformation of abstract spin networks. 

6 Outlook 

In this review we discussed in detail non-Abelian Chern-Simons theory for arbitrary gauge 
group on a three-dimensional manifold with boundary. We examined the relation of this field 
theory to the WZW conformal field theory in two dimensions and saw that, after canonical 
quantization, the wavefunctions of the Chern-Simons theory are the conformal blocks of the 
WZW current algebra at level k. By using the Chern-Simons action as the measure of a path 
integral which defines the expectation value of the Wilson loop, we described Witten's intrin- 
sic three-dimensional definition for the Jones polynomial. We then described a reformulation 
of general relativity on an 5*0(3) Yang-Mills phase space, and defined the fundamental loop 
representation of the non-canonical T-algebra for A = which led to a purely combinatorial 
structure called the spin network. For A^Owe solved the constraints to obtain an exact state 
of quantum gravity whose classical limit is de Sitter spacetime. The presence of A led to the 
framing of loops and therefore to a quantum deformation of the spin networks, resulting in a 
discrete spacetime of punctured Ricmann surfaces with finite genus. 

In loop quantum gravity, we have understood that via the loop transform one can represent 
states of quantum gravity via Wilson loops (and integrals of Wilson loops over the underlying 
gauge field A), and hence by the geometry of knots and links embedded in the three-space. 
The loop transform associates a functional integral to a given loop in three-space, and it is seen 
that in the case of the Chern-Simons state that this transform is, up to framing, a topological 
invariant of the loop (as discussed in more detail below) . One knows that these same topological 
invariants arise by taking bracket evaluations of the knots at appropriate values of the bracket 
variable A. As a result, it is a well-grounded move to shift from the functional integral approach 
to an approach that evaluates knotted and embedded spin networks in three-dimensional space. 
Each such evaluation can be regarded as a specific evaluation of a quantum gravity loop trans- 
form, and one can look for geometric and physical intepretations of these evaluations. One can 
think of each loop state as having all its geometry concentrated along the loop in the fashion 
of a three-dimensional embedded delta function. The metric that arises here assigns an area 
measurement to any test two-dimensional surface that encounters the spin network loop. Each 
intersection point gives the value Stt^HGnC 3 y/ j(j + 1), where j labels the spin network edge, 
and we have retained the speed of light c and Planck's constant ft. One adds up all the contri- 
butions to the area coming from the different spin network lines in a given or chosen family of 
embedded spin networks in three-dimensional space. The remarkable point about this method 
is that in the context of loop quantum gravity, the area measurements are quantized in this 
fashion. The fundamental quantum of area is given by a surface with a single puncture, with 
the edge carrying a spin j = 1/2. In this case A = 4V3irjhG N c- 3 w 1.808 x 10" 69 m 2 . Surfaces 
make quantum jumps between states in the spectrum of the area operator. This, in turn reflects 
on many other issues such as the discrete nature of space and time at the Planck scale and the 
information content of black holes. For more information on these and related matters we refer 
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the reader to [43-46,61-65]. 

Having discussed the achievements of the approach to quantum gravity presented in this 
review, we will now turn to a brief discussion of several key problems that are still open. These 
include: 

• The physical Hilbert space: Time evolution is governed by the Wheeler-de Witt equa- 
tion Sip = 0, where S is the operator corresponding to the complicated Hamiltonian 
constraint in (85). The theory should express four-diffeomorphism invariance but the 
canonical decomposition requires a choice of external time parameter, thus breaking gen- 
eral covariance. The states described here form the "kinematical" Hilbert space; that is, 
these states only solve the quantum constraints Q 1 and Vj. What is needed is a set of 
functional that will solve the Wheeler-de Witt equation and define the physical Hilbert 
space, and on it a physical inner product. A possible candidate for the Wheeler-de Witt 
operator has been defined by Thiemann in [66], but it is not obvious that this operator 
in the classical limit describes general relativity at all. 

• Reconstruction problem: Any candidate theory of quantum gravity should in some limit 
reproduce a classical background metric that satisfies the Einstein equations. Further- 
more, we expect that a good theory of quantum gravity coupled to matter should in some 
low-energy limit reproduce quantum field theory in curved spacetime. This is, however, 
a very non-trivial problem to address. Matter has been coupled to gravity via the self- 
dual variables and has been quantized [67-71]. The problem has only been solved on the 
kinematical Hilbert space. As mentioned above, an inner product for physical states is 
currently lacking in loop quantum gravity. Therefore the dynamics of matter coupled to 
gravity, as for example transition amplitudes to describe scattering processes within this 
framework is currently still out of reach. 

• Normalization of the Kodama wavefunction: An analogous wavefunction to the Kodama 
state exists for SU (2) Yang-Mills theory given by 



where g is the coupling constant and Yes is given in (117). This state is a solution to 
[Ef-iB?]V(A) = with Ef = ig 2 S/SA i a . Witten [72] pointed out that the physical norm 
of this theory, which is given by 



is infinite because the configuration space is unbounded from below. He conjectured 
that this should be true also for the Kodama state. This issue is addressed in [73], 
where it is argued that this argument does not extend to gravity because the Kodama 
state, as any other states that are physical states, are expected to be nonnormalizable 
in the kinematical inner product. The Yang-Mills theory above is not diffeomorphism- 
invariant and hence the state (124) is nonnormalizable in the physical inner product. In 
[73], in order to address the issue of normalizability in the physical inner product, the 




(124) 




(125) 
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authors studied the linearized Kodama state over a de Sitter background. They found 
that although this linearized state is delta-function normalizable in the Euclidean theory, 
it is not normalizable in the Lorentzian theory. This raises an important issue as to 
whether the full Kodama state discussed in this paper is a physically reasonable state to 
describe the ground state of quantum gravity with a positive cosmological constant. 

Most of these issues may be solved in the future using the state sum and spin foam models 
for quantum gravity that are currently under intensive investigation. These are based on a 
covariant path integral approach and hence do not break general covariance. There is not room 
in this short review to do more than mention these aspects. We refer the reader to [74-83] for 
more information on state sum models and [84-88] for more information on spin foam models. 
Several important problems have now been addressed in this direction: 

• The reconstruction problem has begun in [89, 90], employing renormalization group meth- 
ods of statistical physics to find the low-energy limit of spin foam models. 

• Work is in progress on spin foam models of matter coupled to quantum gravity [91-93]; the 
question of particle scattering within the framework of background-independent quantum 
gravity has been addressed in [94]; an attempt has been made at obtaining the graviton 
propagator via spin foams in [95] . 

• Recently the idea has been put forward, that perturbative quantum gravity can be treated 
as an expansion around topological BF theory [96] . This is related to the fact that gravity 
can be written as BF theory with the constraint that the B field comes from the metric. 
An implementation of this idea has been applied to standard Yang-Mills theory in flat 
spacetime [97] and was shown to be a viable technique for perturbation theory of generally 
covariant theories. 

We conclude with a few remarks on the background dependence of matter quantum field 
theories. There it is always assumed from the outset that the background spacetime is con- 
tinuous to arbitrarily small scales. In the non-perturbative quantization of pure gravity we 
generically end up with a discrete space with length, area and volume eigenvalues. So the ini- 
tial assumption of quantum field theory is false in this context, and may in fact be the origin of 
the ultraviolet divergences that plague the quantum field theories. Indeed in his QSD V article, 
Thiemann [98] shows that the divergences can be fully accounted for if we neglect the discrete- 
ness of space, and has argued that the Planck scale may provide a natural cutoff regularization 
of the Hamiltonian constraint of matter quantum field theories. We can go further. By looking 
at the modifications to matter interactions after taking the discreteness of space into account, 
we may find a better understanding of nature, built on a well-defined set of principles. 
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A Conventions 

Here we collect the conventions used in the paper. Unless otherwise stated, we work on a 
four-dimensional Lorcntzian manifold M. which admits constant time hypcrsurfaces S. We use 
the following index conventions: lower-case greek letters \i,v,... € {0,1,2,3} are spacetime 
indices; lower-case roman letters € {1, 2, 3} are spatial indices; upper-case roman letters 

I, J,... G {0, 1, 2, 3} are internal SO (3) tangent space indices; upper-case roman letters with 
hats I, J, . . . e {1, 2, 3} are gauge-fixed internal SO (3) tangent space indices. Spacetime indices 
are raised and lowered with the spacetime metric g^ vi spatial indices are raised and lowered 
with the induced spatial metric q t j on a constant time slice, and internal indices are raised and 
lowered with the Minkowski metric r]u = diag(— 1, +1, +1, +1). Here spacetime and spatial 
indices are consistently written before the tangent space indices for all quantities on A4. The 
Einstein summation convention is used for repeated upper and lower indices. Units are employed 
for which the speed of light c and Planck's constant h are unity, and Newton's constant Gn is 
retained. 

B Hamiltonian Formulation of Classical Mechanics 

For convenience we summarize here the Hamiltonian formulation of a classical unconstrained 
system with a finite number of degrees of freedom [99,100]. We consider a system that is 
governed by the action S[q,q;t] = J L(q,q;t)dt, where t is an integration parameter with an 
overdot denoting d/dt, and L is known as the system's Lagrangian. This is a functional, and 
is defined as the kinetic energy minus the potential energy. In most cases the kinetic energy 
T = T(q) is a function of the system's velocities q = q{t) = {<z Q (i)}a =1 and the potential 
energy V — V(q) is a function of the system's coordinates q = q(t) = {<7 a (£)}„ = i- The set of 
all coordinates {q a }a=i are said to define an n dimensional manifold called the configuration 
space Q. However, the defining action S[q, q; t] exludes the more interesting gauge theories 
because the Lagrangian only describes systems with a finite number of degrees of freedom. 
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For the former we would need to consider n — > oo. Furthermore, the parameter t in most 
nonrelativistic systems is the time which measures evolution, but this breaks down when we 
consider relativistic motion since time is no longer a parameter. 

Now, the functional S = S[q, q; t] is diffcrentiable, and the equations of motion of the system 
are obtained by demanding that the action be stationary with respect to variations of q and q. 
We find that 



^dq a dq a J J \dq a dt dq a , 

Here we integrated by parts using 5q a = d(Sq a )/dt. The boundary term vanishes because we 
are taking the variation with respect to fixed endpoints. If we impose SS = for arbitrary 5q a , 
then it follows that 

dL d dL n 

Oq^-dtW= - (12?) 

This is a set of n second order differential equations in q a , q a , and t. These are called the Euler- 
Lagrange equations of motion, and by Hamilton's principle of least action (which states that the 
motion of a mechanical system coincides with the extremals of the functional S[q, q; t]) determine 
the (time) evolution of the system. Furthermore, the canonical (generalized) momenta of the 
system are given by 

dL , N 

Pa = W' 1 (128) 

these are of central importance in the Hamiltonian formulation to be introduced below. They 
are generalized momenta in the sense that they represent a generalization of the Newtonian 
notion of momenta given by p = mq. Equations (127) can be shown to reproduce Newton's 
second law F = ma in the form dp/dt = — VV. 

An equivalent approach to Lagrangian mechanics outlined above is to use a Legendre trans- 
form to obtain a set of first order differential equations to replace (127). To this end, we let 
L = L(q,q;t) and p a — dL/dq a . Now, we define the Hamiltonian function H = H(q,p;t) 
through the Legendre transform 

H(q,p-t) = q a p a -L{q,q;t). (129) 

The differential of this function is 

dH = q«dp a +p a dq a -^dq«-^dq a -^dt 

dL 

= q a dp a -p a dq a - —dt. (130) 

In general, however, we should expect the Hamiltonian H — H(q,p;t) to have differential of 
the form 

fszj f)FT f)M 

"^W^^dp-^^- 
Comparing this to (130) we obtain the Hamilton equations of motion 

pi TT Ci TT 

q a = jr- and Pa = -° (132) 
dp a dq a 
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for which the Lagrangian must satisfy dL/dt = —dH/dt. If the transformations connecting the 
rectangular and canonical coordinates are independent of time so that the kinetic energy is a 
homogeneous quadratic function of the q a , and the potential energy is independent of velocity 
then it follows that the Hamiltonian is equal to the total energy of the system and is a conserved 
quantity provided it does not explicitly depend on time. In fact, the Hamiltonian is a conserved 
quantity provided that the system is closed, i.e. a system which docs not interact with its 
surroundings. 

Let's now consider a function denoted by A = A(q,p;t), for which we calculate the total 
time derivative to be 



dA ( dA dq a dA d Pa \ dA 



dt \dq a dt dp a dt ) dt 

dA dH dA dH\ dA 



dq a dp a dp a dq a J dt ' 



(133) 



Here we made use of the Hamilton equations (132). This can be written more abstractly if we 
define the Poisson bracket 1 



dq a dp a dp a dq a 



so that 



dA , , dA , . 

M= {A > H} + lH- (135) 

In most cases A will contain time dependence implicitly so that dA/dt = 0, and thus we have 

dA 

= {AH}. (136) 

This beautiful expression contains in it all of classical mechanics. If we evaluate the Poisson 
bracket for q a and p a we obtain the equations (132). In this form, we see that the Hamiltonian 
is the generator of time translations. The set of all coordinates {q a , p a }a=i are sa id to define 
a 2n dimensional manifold called the phase space P. A point (q a ,p a ) on T is called a state, 
and a function F = F(q,p;t) on T that satisfies Hamilton's equation F = {F,H} is called an 
observable. The set of n second order equations on Q have been replaced by a set of 2n first 
order equations on T. 

A more absract formulation can be obtained if we define the phase space as a dual tangent 
bundle, or cotangent bundle, with coordinates (q a ,p a ) on which the Hamiltonian H : T*{N) — > 
R is defined. This can be seen in the following way. Recall that the Lagrangian L = L(q, q; t) 
is defined over an n-dimensional manifold called the configuration space Q. The q a are the 
components of a vector at the point (q 1 , . . . ,q n ). Thus the Lagrangian L : T(N) — > R is a 

1 The Poisson bracket satisfies the following identities: 

1. {F, G} = —{G, F] 

2. {F + H,G} = {F,G} + {H,G} 

3. {FH,G} = F{H,G} + {F,G}H 

4. {F, {G, H}} + {G, {H, F}} + {H, {F, G}} = 
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function on the tangent bundle. Now the Hamiltonian is obtained from the Legendre transform 
(129). The first term q a p a can be visualized as a pairing of an element of the tangent space 
with its dual. In replacing q a with p a as the new set of independent variables, the Legendre 
transform replaces the tangent bundle T(N) with the cotangent bundle T*(N) and thus defines 
the Hamiltonian as a function H : T*(N) — ► R. This explains why in the Einstein summation 
q a p a = J2 a QaPa for some canonical set of variables {q a ,p a }2=i we write the position q a with 
contravariant index and the momentum p a with covariant index. 

The connection between conserved quantities and the symmetries of a Lagrangian is con- 
tained in Noether's theorem. The essence of this theorem is that, if the Lagrangian is unaffected 
by a transformation that alters one of its coordinates, then the Lagrangian is invariant, or sym- 
metric, under the given transformation, and the corresponding current is then conserved. In 
particular, there are seven constants associated with the motion of a closed system: homogene- 
ity of time implies that the Lagrangian is invariant under time translations, which leads to the 
conservation of energy; homogeneity of space implies that the Lagrangian is invariant under 
spatial translations, which leads to the conservation of linear momentum; isotropy of space im- 
plies that the Lagrangian is invariant under spatial rotations, which leads to the conservation 
of angular momentum. In the Hamiltonian formulation, a conserved quantity is simply one for 
which the Poisson bracket vanishes, i.e. dA/dt = {A, H} = 0. This provides a simple proof of 
the conservation of energy since {H, H} = must be satisfied by the antisymmetry property of 
the Poisson bracket. Finally, we define a cyclic coordinate q a to be one that does not appear 
explicitly in the Lagrangian. Then (127) reduces to p a = 0, in which case the Hamiltonian 
does not explicitly depend on q a either. Thus the canonical momentum conjugate to a cyclic 
coordinate is conserved, i.e. {p a ,H} = 0. 

C ADM Decomposition 

In order to cast the Einstein-Hilbert action S[g] = (1/16ttGn) f d 4 Xy/g(R — 2A) with R = R[g] 
into Hamiltonian form, we need a time with which we can define the conjugate momenta. 
To this end, we split the spacetime into a space and a time. We assume that M has the 
topology [0, 1] x S where S is an open or closed three-dimensional surface; this is a segment of a 
universe between initial surface {0} x £ and final surface {1} x £ both assumed to be spacelike. 
We also assume the existence of constant time hypersurfaces {St} endowed with coordinate 
systems {x 1 } and induced metrics {qij(t, x )}. After an infinitesimal time translation from S t 
to St+dt, the change in proper time will be dr = Ndt, and this leads in general to a spatial 
shift x\t + dt) = x l (t) - N l dt. Here N = N(t,x l ) is the lapse function and N l = N^t.x 1 ) is 
the shift vector. The spacetime interval between two points (t, x l ) and (t + dt, x l + dx l ) on M. 
is then 

ds 2 = -N 2 dt 2 + qij (dx l + N'dtXdx* + N 3 dt). (137) 
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Figure 4: The line element expressed in the ADM decomposition. 
This is shown in Fig. 4. In components this gives the metric tensor and its inverse 

9v* = AT ( 138 ) 

(_ 1 N 3 \ 

Here q lJ is the inverse of qij. From the form of the inverse metric above we note that q lJ are 
not the spatial components of g^ u . 

Using the induced metric we can determine the curvature tensor R Z j kl that is intrinsic to the 
hypersurfaces St. This is determined by parallel transporting tangent vectors on St. However, 
we can also define an extrinsic curvature to describe how the hypersurfaces curve with respect 
to the embedding spacetime. This is described by the behaviour of vectors normal to S t . For 
n a = (N, 0, 0, 0) a unit normal to S t in the lapse and shift decomposition above, this is given 

by 

Kij = - (3) ViNi VjNi), (140) 

and with this we find that the Einstein-Hubert action becomes 

S = — ^- [dt [ d 3 x^qN(WR -2A + -(qijKV) 2 ) 
LbirGN J Js 

+boundary terms. (141) 
Details can be found in [101, 102]. The canonical momenta can now be read off from the action: 

1^- = ^q(K^ - q^(q kl K kl )) (142) 
oqij 



tt = ^=0 (143) 
SN 

tt 1 = ^- = 0. (144) 
SNi 
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There are four primary constraints tt° and n l . The first momentum can be inverted to give the 
extrinsic curvature 

= ^-{-K ii -q ij (q k m kl )) (145) 

and hence the action (141) becomes 

S = JdtJ d 3 x{ir ij qij — NH - N,H l ) (146) 

H = -^=(7r«7r« - (qijir ij ) 2 ) - V?( (3) i? - 2A) « (147) 
H i = _ 2 (3) Vj7r u „ o. (148) 

Here, the Hamiltonian constraint H and momentum constraints H l are secondary constraints, 
and the lapse and shift are Lagrange multipliers. The action (148) is in the form of a Lagandre 
transform / p^q^ — H for which the term 7r 4J ' qij plays the role of the term Thus we 

identify the quantity Hadm — NH + NiW as the gravitational Hamiltonian (called the ADM 
Hamiltonian for Arnowitt-Deser-Misner) . Because the Hamiltonian is a linear combination of 
constraints it vanishes on-shell (i.e. when the equations of motion are satisfied). From Hadm 
we find the equal-time Poisson brackets 

{ qij (x),<ir kl (y)} = torG N $6 l j + 5\6!?)6 3 (x - y) (149) 
{qij(x),q k l(y)} - (150) 
{7r«0r),7r fcI (i/)} = 0, (151) 

and / d 3 x5 3 (x — y)f(y) — f(x) for any function f(x). 
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